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Abstract. Logarithmic perturbation theory (LPT) is developed and appHed to 
quasinormal modes (QNMs) in open systems. QNMs often do not form a complete set, 
so LPT is especially convenient because summation over a complete set of unperturbed 
states is not required. Attention is paid to potentials with exponential tails, and 
the example of a Poschl- Teller potential is briefly discussed. A numerical method is 
developed that handles the exponentially large wavefunctions which appear in dealing 
with QNMs. 
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1. Introduction 

1.1. Logarithmic perturbation theory 
Eigenvalue problems of the type 

H(j) = X(j) (1.1) 

occur in many branches of physics; here A may be related to the frequency uj by 
X = u (Schrodinger equation), or by A = cj^ (wave equation or Klein-Gordon equation). 
Perturbation theory is useful for systems that depart slightly from an ideal solvable 
configuration. Apart from the standard Rayleigh-Schrodinger perturbation theory 
(RSPT), a useful alternative focuses not on itself, but on its logarithmic derivative 
/ = 4>'/(f). Known as logarithmic perturbation theory (LPT) |1|, ^, 0, ^, this method 
is commonly applied to 1-d bound state problems, especially the ground state. For 
excited states, one either has to first factor out the zeros [Q, or detour around them in 
the complex plane LPT has also been developed for bound states in 3 d LPT 
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avoids sums over intermediate states, and comparison with RSPT can lead to useful 
sum rules 0. 

The bound states, or normal modes (NMs), are solutions with — ^ as x ±oo. 
Other boundary conditions are also important physically. Scattering states and the 
phase shift can be handled by LPT using only on-shell information ||^. In this paper, 
we develop LPT for wavefunctions that are outgoing at infinity — quasinormal modes 
(QNMs). 



1.2. Quasinormal modes 

In conservative systems, NMs are factorized solutions = e~*'^*(/)(a;) with cf) 

satisfying an eigenvalue equation such as and nodal boundary conditions at 



X — > ±oo. The counterparts in open systems are QNMs; these factorized solutions 
satisfy outgoing wave boundary conditions at s ^ ±oo, so that Imtu = —7 < 0. 

QNMs are important from many points of view. A laser is often discussed in 
terms of its "modes", i.e., the spectral lines with finite widths 7, which are precisely 



these QNMs [|T0[- Quantum-mechanical resonances are likewise central to scattering 
111 , |12|, 0, |14[, and as intermediate states in high-order transitions. Gravitational 
waves from the vicinity of a black hole are likely to be detected in the next decade 
by facilities such as LIGO and VIRGO |]T5|. The radial wavefunction describing the 



propagation of gravitational waves in any angular momentum sector satisfies a Klein- 
Gordon equation with a potential V{x) [|T^, |T^. Theoretical studies [|T^, |T^ show that. 



at least for an intermediate time domain, the waves are dominated by a ringing signal, 
which is readily identified as the superposition of QNMs []TB|, |2D[ . If the relationship 
between the characteristics of the ringing signal (i.e., the QNMs) and the spacetime 
curvature could be better understood, gravitational waves have the prospect of becoming 
a novel astronomical probe. In these cases, the background is a Schwarzschild metric 
plus perturbations (e.g., due to an accretion disk), so perturbative treatments will be 
useful. 

In the present context, three properties of QNMs should be emphasized. First, 
their numerical determination is notoriously difficult. This is most simply seen in the 
"shooting" algorithm: choose an arbitrary cj, integrate from one end (say x = for a 
half-line problem or a full-line problem with definite parity), identify the coefficient of the 
"wrong" solution at the other end (say x 00), and vary uj until this coefficient is zero. 
For NMs, the exponentially large "wrong" solution is readily identified. For QNMs, the 
"wrong" solution, which is 0{e~'^'^^) relative to the "right" solution, is difficult to extract, 
especially when 7 is large. The numerical difficulties make perturbation methods even 
more relevant than would be the case for NMs. 

Secondly, RSPT is inapplicable for two reasons. Its usual derivation rehes on the 
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hermiticity of the system, which is now lost. Moreover, because the QNMs are in general 
not complete, one cannot sum over intermediate states. Even in circumstances where 
the QNMs turn out to be complete ]T2|, |T|, g^, g2], |2|, ^ |6[ , a scheme such as LPT 
would still have definite advantages, because it makes no reference to the higher states 
with large 7. 

Thirdly, it is readily shown that any QNM can have at most one node on the real 
X axis. Except for the origin for the odd-parity sector of symmetric potentials, there 
is no reason why any root of the complex equation 0(a;) = should lie on the real x 
axis; those cases that do are therefore "accidental" in the sense that they occur only for 
specific values of the parameters defining the potential V{x) — in other words on a set 
of measure zero in parameter space. Therefore the nodal problem which plagues LPT 
for excited NMs is here absent generally. 



1.3. Outline of paper 

Section ^ develops LPT for QNMs, and discusses the generalized norm that emerges 
as a result. The most explicit general form for the second-order correction, together 
with an illustrative example, are given in Section 0, focusing on those cases where both 
the original potential and the perturbation have finite support. The situation becomes 
slightly more complicated if the potentials have tails, and the case of exponential tails 
is discussed in Section ^. A conclusion is given in Section ^. 



2. Perturbation theory 



2.1. Formalism for the eigenvalue 

We deal with the Klein-Gordon equation: 

'dl - V{x) + uj^] (f){x) = . 

The Schrodinger equation is included by simply re-labeling uj^ 
derivative f{x) = (j)' (x) / (f){x) satisfies the Riccati equation 



f'{x) + f{x)-V{x)+uj' 



. 



(2.1) 

uj. The logarithmic 
(2.2) 



We let / denote the logarithmic derivative corresponding to an eigenvalue, so that 
it satisfies the two boundary conditions f{x) ±iu as x — > ±00. At a general 
frequency, however, we can define similar functions (j)±{uj,x) and their logarithmic 



derivatives f±{uj,x) as solutions to (2J.) and ( p.2|) , but with each function satisfying 
only one boundary condition, namely f±{uj,x) — > ±iuj as x — > ±00. At an eigenvalue, 
/- = /+ = /• 
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Now let the potential be perturbed 

V{x) = V^{x)+^iVi{x), (2.3) 

where /z is a formal small parameter. The eigenvalue uj and the function / are both 
written in powers of /i f: 

UJ = uoq + ^uoi + H , (2.4) 

/ = /o + /i/i+/^V2 + --- = /o + ^7, (2.5) 



where /o, assumed known, satisfies the Riccati equation (2^) with the potential Vq and 
frequency ujq. 

Now divide the real line into three regions (— oo,L_), (L_,L+) and (L+,oo). If 
the original potential and its perturbation both have finite support within the central 
interval, then the asymptotic regions are trivial, and the simplest examples will be of 
this type. 

First consider the central region, and put ( |2.4| ) and ( |2.5| ) into ( ^.2| ). Upon comparing 
powers of /i, one finds 

fn + 2/o/n + 2^o^„ = K , (2.6) 



for n = 1,2, in which Vi is the perturbing potential in ( p.3| ), and V^, > 1, is 
a shorthand for the following combination of lower-order quantities, to be called the 
effective nth order potential 



n-l 



Vnix) = - XI [fi{^)fn~i{x) + UJiUJn-i] ■ (2.7) 



i=l 



Using the integrating factor exp[2 dy fo{y)] = (l>l{x), one gets from (|2.6|) 



/„(x)0o(x) = / dx [Vn{x) - 2u;ou;„] 4>l{x) . {2.t 



We now need to match the central solution to the two asymptotic regions. Assume 
that the latter have been solved with outgoing wave boundary conditions at spatial 
infinity, and denote the logarithmic derivatives to be matched as 

D±{uj) = f±{uj,L±). (2.9) 

t It is a property of LPT that one need focus only on one state at a time. Therefore, a label for 
different QNMs will in general be suppressed. 
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Note that i5± will contain two types of changes from the unperturbed case. First, at 
fixed cj, the wavefunction when integrated inwards from ±00 will suffer changes because 
of Vi{x) in the two asymptotic regions; these are expressed through 

D±{u) = D^oiuj) + fiD^iH + fi^D±2iuj) + ■■■ . (2.10) 

Secondly, there will be changes because the value of u itself shifts according to (|2.4|) . 
In particular, the exact logarithmic derivative is f±{u!,L±) = D±{uj), whereas the 
corresponding unperturbed quantity is fo{L±) = D±q{ujo). Thus 

9{L±) = f±{uj, L±) - /o(L±) = D±{uj) - /^±oM ■ (2.11) 

Using ( |2.4|) and (|2.1CI|) and developing the right hand side of ( p.ll| ) in powers of /i, we 
can find that /„ in ( p.8|) should be matched to 

/„(L±) = c^nZ^^oK) + A±„ , (2.12) 

where A-|-„ does not contain ojn] exphcitly but in shorthand 



Ai = A , 

^2 = D2 + ujiD[ + ]^ujIDI, (2.13) 

etc. In the above, the subscripts ± have been omitted from all quantities, and all Dn on 
the right are to be evaluated at ojq. In short, one requires a knowledge of the perturbation 
in the asymptotic region (D„, n > 0), as well as a knowledge of the unperturbed problem 
slightly away from the original frequency (derivatives of Dq). 

Putting these into ( p.8|) and collecting terms involving one finds the central 
result 

_ (0o|K|</'o) , . 

^•^ ~ o — 7T7TT ' \^-^^) 
2t^o(0o|0o) 

in which we have introduced the suggestive notation 



(0o|K|0o) 
(0o|0o) 



rfxV;(x)0o(x) + -A+„0^(L^ 



A_„0g(L. 



i+ 1 



2u;o 



(2.15) 
(2.16) 



This expresses the nth order correction in quadrature in terms of lower-order quantities 
(provided the asymptotic regions have been solved to give /S.±n and D'^q). 
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The division into three regions is arbitrary, and the whole expression ( p.l4[ ) must 
be independent of L±. Moreover, the numerator and denominator must be separately 
independent of L±, because the numerator depends on the perturbation, whereas the 
denominator relies only on the unperturbed system; an explicit proof can be constructed 
by calculating d{(l)o\(f)o) /dL+, and then using (^.21). 

Thus, in both (2.15) and (2.16), we can formally take L± ±oo and write 

/oo 
dxVn{x)^lix) , (2.17) 
-oo 

/oo 
dx^lix) . (2.18) 
-oo 

These formal expressions do not converge; (2.15) and (2.16) may be regarded as ways of 
regularizing them. In Section ^ we shall discuss various different ways of giving meanings 
to these formal integrals. 

Evidently, the numerator should be regarded as a generalized matrix element, and 
the denominator should be regarded as a generalized norm. We now develop this 
interpretation. The corrections to the eigenfunction will be given in Section |2.3| . 



2.2. Generalized norm and matrix element 

For the simplest case where Vo{x) vanishes outside the interval {L^,L+), the solutions 
in the two asymptotic regions are exactly e^*'^^, and D±q{uj) = ±iu!. The generalized 
norm (2.16) simplifies to 



(0o|0o) 



dxi 



(x) + 



(2.19) 



In this form applicable to potentials without tails, the generalized norm has been 



introduced previously both for the wave equation [0, the Schrodinger equation 
28| and the Klein-Gordon equation p4|, and its properties discussed. It has been 



shown to be equivalent to another form first given by Zeldovich p9[ , which did not have 
the surface terms, but instead required a process of regularization |[Tl|, |l^, 13 1 which 
is less convenient for actual evaluation (especially numerical evaluation). The present 
result, in the more general form (2.16), is however applicable to potentials with tails, 
and examples will be given in Section ^. 

We next briefly describe the properties of this generalized norm, and argue why it 
deserves to be so named. 

First of all, suppose the system parameters can be tuned so that the leakage of 
the wavefunction approaches zero (e.g., if Vq{x) contains a tall barrier on both sides). 
Then L± can be chosen so that 0o(-^±) ~ 0; the expression in (2.16) then contains 
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only the integral. Moreover, when the leakage is zero, the frequency is real, and the 
wavefunction has a constant phase, which can be chosen to be real; thus 0q = |0oP- 
The expression (2.16) then reduces to the usual (real and positive-definite) norm for a 
NM. Because of this property, and because it appears in the denominator in ( p.l4| ) to 



scale the wavefunction, it is appropriate to call this quantity the generalized norm. 

Nevertheless, it has some unusual properties, (a) It involves 0q rather than |0oP, 
and is in general a complex quantity, (b) It involves a surface term, though the value 
of the entire expression is independent of the choice of L±. Thus, it is not a genuine 
norm, and the term is merely a shorthand for "a bilinear map that appears in the place 
of the norm in perturbation formulas such as (|2.14|) " . 

It is hardly surprising that perturbative results are expressed in the form of a matrix 
element divided by a normalizing factor, as in (p.l4|) , but it would not have been obvious 



what the normalizing factor should be. The point is that for a QNM, the wavefunction 
behaves as cI)q{x) ~ e*"^"^ oc e"'^^ as x — > cxo and 70 = — Imcug, so that an expression such 
as (2.18) (and even more so for the analogous formula with |(/>oP) would be divergent. 
An important result of the present paper is that we give a precise way of normalizing 
such QNM wavef unctions. 

We have already remarked that the generalized norm is not real, and neither is the 
(diagonal) matrix element. Far from being a problem, this is necessary, in that the 
result ( p.l4|) gives the corrections to both the real part and the imaginary part of the 
frequency. Thus, despite the formal similarity to the analogous problem for NMs, the 
present formalism in fact contains twice the amount of information. 

Some of the properties above, in particular the validity without regularization of the 
simpler form (2.17), relies on Vn behaving mildly at infinity. It is therefore appropriate to 
demonstrate that if the perturbation Vi has finite support, then so does all the effective 
potentials generated via (^3)- Consider for simplicity only the asymptotic interval 
(L_|_,oo). Now the exact logarithmic derivative is iuj, whereas the unperturbed analog 
is iojQ. This then gives 

g^ = -{uj-ujof = -[Y,li'uA . (2.20) 



The nth order term in the above expression then ensures that Vn in (|2.7| ) exactly vanishes 
in this region. 



2. 3. Wavefunction 



To complete the iterative procedure, we also need the eigenfunctions. This can be readily 
obtained by integrating ( |2.8| ) to an arbitrary point, and using ( |2.11| ) as the boundary 
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condition, which yields 



fn{x)(j)l{x) = [uJnDioiuJo) + A_„J (j)l{L.) + / dy [Vn{y) - 2uJoUJn] (flijj) • (2.21) 



One could write an alternate expression using the boundary condition at Lj^ and 
integrating from the right. Consistency is guaranteed if uJn has been correctly evaluated 

by (EH- 

Thus we have in principle an order-by-order iteration scheme for the QNMs; namely, 
use (p.l4|) to obtain cji, and (|2.21|) to get /i; this is then put into (|2.7|) to find V2, etc. 



3. Explicit form of higher-order corrections and an example 
3.1. Higher-order corrections 

The perturbative formulas would be more useful if they could be written explicitly rather 
than recursively. In general, the nth order correction to the frequency must take the 
form of an integral over V{xi) ■ ■ ■V{xn)\ moreover, the perturbing potential can only 
act if it is "sampled" by the wavefunction 0o(^)- ^^^^ ^^^^ turn out to be convenient 
to remove a constant from Vi, and we are led to define 



W{x) = [V,{x)-2uo^uo,]<Pl{x) . (3.1) 



The constant subtracted renders the integral of W zero (see ( p.l4|) ) 



Thus we expect to be able to write the nth order correction in the form 
If " 

^n = - 7— 7— T dXi- ■ -dXnS \\W{Xj)'^n{Xl,- ■ ■ ,Xn) , (3.2) 

2t^o(0o|0o) J jJi 

where 5* = 0{xi —X2) ■ ■ ■ 6{xn-i —Xn) makes use of the symmetry among the coordinates 
Xi, - ■ ■ ,Xn to restrict the integration to one sector {9 is the unit step function), and the 
weight function constructed out of (po, scales as (0q)^~". 

We now try to write out 002 in essentially this form, and evaluate the weight function 
"^2- For this purpose we consider the simpler case of a half line < x < 00, with the 
potentials satisfying Vq{x) = Vi{x) = for x > a, and all 0(x = 0) = 0. (This may 
be regarded as the odd-parity sector of a symmetric problem.) Thus, all the surface 
contributions at L_ are eliminated, while at the right hand side we may take = a 
and D'_^_Q = i. Thus, the generalized norm is 



a 



dx(f)Q{x) + — (/>o(a) ■ (3.3) 
2uJo 



By using i \2.1\ ) for V2, we can write the second-order matrix element as | 



(0o|V2|0o) = - / dx(f)Q^{x) /i(x)0o(x) -ujf dx(f)l{x) . (3.4) 



Now from ( |2.21| ) we have /i(x)0g(x) = Jq dyW{y). Putting this into (^.41), changing 



the order of integration and also using ( p. 3D to simplify the second integral in (p.4|) then 
leads to 

2u;o(0o|0o) 2a;o 4u;o (0o|0o) 

in which the weight function is found to be 

^2{y,z) = ^2{y) = -2 rdx<pf{x) . (3.6) 

Jv 



Thus, except for the last two terms in ( |3.5| ), which do not involve an integral, the second- 
order correction has been cast in the form (p.2|) , which is the most explicit form possible 
for an arbitrary perturbation. 

3. 2. Example 

We now illustrate these formulas by a very simple example. Let the unperturbed system 
be defined by a step 

VQ{x) = VQe{b-x) , b<a. (3.7) 

The unperturbed eigenfunctions are 

^ r ^ _ / ^singx , x <b , 

\ Asing6e-o(--^) , x>b, ^ ' 

where the condition of outgoing waves results in the eigenvalue equation for q: 



q cotqb = iJ q^ + Vq = iujq . (3.9) 



Here loq is the unperturbed frequency. There are of course many solutions to (|3.9| ), and 
we pay attention to any one of these. 

The norm is readily evaluated by (|3.3|) to be 



A% ( sm2qb sin^ qb tanqb\ 
(0o|0o) = — (^1 - ^b ) ■ ('-'^^ 

I For QNMs, the wavefunction does not have any nodes apart from the one at a; = imposed by the 
boundary condition. 
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For any perturbation Vi with support on (0, a), the first-order shift is then 

/o dx sin^qx Vi (x) + dx sin^gfo e^^-^o^^-*) Vi (x) 

^1 — 7 \ ) (3-11) 

ujob (l — sm2qb/2qb — sin qb tanqb/qb) 



while the second-order correction is given by ( p.5| ), with the weight function being 

^ I {2/A^){cotqb - cotqy) + C [l - e2-o(^-)] , y<b, 

•^yy) I ^(g-2ia.oy_e-2*c.oa^ , y>b, ^' ^ 

where C = i/{uoQA^ sm^qb). 

These formulas then allow the corrections for any perturbation Vi to be obtained 
by direct quadrature; more importantly, they exhibit how the perturbation Vi acts to 
shift the complex eigenvalues. 

To be specific, let the perturbation be a bump of width w centered at a position xq: 

{w w 
, otherwise . 

Figure 1 shows the trajectory of the lowest eigenvalue, u, in the complex frequency 
plane as xq is changed, for fixed Vq = 100, b = 1, fi = 10, w = 0.1 and (a) Vi lying 
inside the interval (0, b), (b) Vi lying outside the interval (0, 6); the exact results (circles), 
first-order perturbation computed by ( |3.11| ) (dashed line) and second-order perturbation 
computed by ( ^.5] ) (solid line) are shown together for comparison. Figure 2 shows the 
magnitude of the remaining error for the unperturbed eigenvalue (solid line), first-order 
perturbation (dashed line) and second-order perturbation (long broken line) versus /i, 
for fixed Vq = 100, b = 1, w = 1 and (a) Xq = 0.3, (b) xq = 1.4. 

Although this example is extremely simple, it illustrates several interesting features. 
First, the remaining error of the nth-order perturbation scales expected. 
This is the case even for the perturbation lying outside the interval (0,6) (Figure 2b), 
in which case the spectrum of QNMs is not complete, and it is not possible to write the 
second-order correction as a sum over intermediate QNMs; this result for the second- 
order correction is testimony to the utility of LPT (as opposed to RSPT). 

Comparison between the two cases in Figure 2a and 2b also reveals that higher-order 
corrections are more significant when the perturbation acts at a more distant position, 
where |0oP is large. 

The behavior in Figure 1 is even more interesting, showing a spiral structure as 
the position xq is changed. Although the perturbation is real and positive, the shift 
can have any phase depending on where the perturbation acts — a situation totally 
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different from NMs in conservative systems. This behavior is most readily understood 
in first-order perturbation theory, for which the general result (|2.14 ) can be written as 



5ui 



5V{x) 



H(x) = 



(0o|0o) 



(3.14) 



where H{x) can be read off from ( p.8|) and (|3.10| ), which makes the spiral structure easy 
to understand. The pattern of the shifts can therefore be very rich, and an example of 
the results for a model astrophysical perturbation of a black hole has been given recently 



4. Potentials with tails 



4.1. The Pdschl-Teller Potential 

When dealing with QNMs, one recurring complexity is the asymptotic behavior 0o(^) ~ 
g27o|a;|^ making norms and matrix elements divergent. As far as LPT is concerned, 
this complication occurs at two different levels. If the potential (and its perturbation) 
vanishes outside a finite domain, then the expression ( |2.19| ) suffices to produce a finite 
expression for the norm, whereas the matrix elements involve integrals only over finite 
domains. Such simplifications also extend to potentials that vanish at infinity faster 
than any exponential. However, when the potential (or its perturbation) decays as an 
exponential or slower (which we shall refer to as a tail), then the evaluation of the norm 
and the matrix element will require more attention. 

In this Section we illustrate the solution of these problems with the example of the 
Poschl-Teller (P-T) potential ^ 

V{x) = Vocosh-\x/b) . (4.1) 



From the point of view of LPT, the P-T potential is interesting because its large 
behavior is exactly exponential: V{x) oc e"^^'/''. The P-T potential, as one of 
a few exactly solvable models, has been studied in depth, in part as a proxy for 
the Regge- Wheeler potential |T^, |TB|, or the Zerilh potential |]T^, which describes 



linearized gravitational waves propagating on a Schwarzschild background. These also 
have exponential tails (as the tortoise coordinate x —oo, i.e., towards the event 
horizon), and consequently their QNMs share certain key properties with those of the 
P-T potential (e.g., a string of QNMs evenly spaced "vertically" in the complex tu-plane, 
— Im co'(j) oc (j -|- ^) §). Therefore a better understanding of the exponential tails may 
be relevant to gravitational waves as well. 



§ The mode index will be indicated in (), but will be suppressed where no confusion arises. 
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The QNM eigenvalues of ([4.1|) are pT 



(4.2) 



where we have assumed 4Vo&^ > 1. The positive (negative) parity sector corresponds to 
even (odd) j. For the purpose of illustrating the LPT formalism, we shall focus on the 
lowest state in each sector, i.e., j = and j = 1. 
Consider perturbations of the width, specifically 



(4.3) 



Because the model is exactly solvable for all b, we immediately obtain the frequencies 
in powers of fi: 



c^o = V^o - i - (j + l)i = a- (j + i) 
1 

4a 



^1 = - - + h)^ ' (4-4) 



etc., where for simplicity we only show the one eigenvalue of each pair with Re > 0. 
We now show how the shift ui can be obtained from LPT. From (14.11) and (14.31), 



Vo(x) = Vo cosh , 

Vi{x) = — 2Vox sinhx cosh~^x . (4.5) 

We show three different ways of handling the divergent integrals (2.17) and (2.18). The 
first two methods are specific to the P-T potential (or other potentials amenable to 
analytic treatment), but these lead to the third method, which is numerical and can 
be applied to any potential with exponential tails. The last method will be the one of 
general interest. 

4-2. Analytic continuation 

Consider for example the j = state. The unperturbed eigenfunction is || 

0o(a;) = (coshx)^'^ , j = 0, (4.6) 

II This expression applies only for uj equal to the eigenvalue; otherwise there is another term with an 
incoming wave, whose coefBcient vanishes at the eigenvalue. The omission of this term does not affect 
the argument based on analytic continuation. 
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where it is understood that uj is to be evaluated at the unperturbed value uj = ujq = 
o — i/2. Take the formal expression (2.18) and define, for any lo for which the integral 
converges, 

/■oo 

N{uj) = / dx4>l{x) , (4.7) 



where 0o is given by ( |4.6| ). The integral N{uj) is well-defined for Im > 0, in which 
domain it is evaluated in terms of the beta function B to be 



N{u) = B{l,~iuj) . (4.8) 

By analytic continuation, this applies to Im u; < as well, and the norm of the j = 
state is then B{^, —1 — ia). Likewise, the wavefunction for the j = I state is 

0o(x) = tanhx (cosha;)*"" , j = 1 • (4.9) 



The same analytic continuation gives the norm as -B(|, — | — icr). 
The matrix element (0o|^i|0o) ioi the j = state is 

-2Vo r dxa;sinhx(coshx)2— 3 = _ v^^o r(l - ^a;) ^ 
j-oo (1 — iuj) r(| — iuj) 

which is convergent even at uj = Uq = a — i/2, and readily evaluated to be 



v/^Vor(l-2a) 



(i-m)r(i-2a) • 

Thus the first-order shift is obtained from ([4.8|) and ( [4. 101) to be 
1 1 
^^ = 4^-2^' 
in agreement with (4.4). 

Similarly, for the j = 1 state, the matrix element is 



-2V, 



dx X sinh^ x ( coshx ) ^ 
Vo 



lUO 



B{l,l-iuj) + 



TT r(l — iuj) 



IUJ 



(4.11) 



where the integral is evaluated for Imcj > 0. Analytic continuation to the eigenvalue is 
required, and gives 



1 

2 L 



B{1 



-la] 



The first-order shift of the j = 1 state is again in agreement with (4.4). 

Analytic continuation, though convenient, applies only when the integrals can be 
evaluated exactly. We therefore present other methods, including numerical evaluation 
of the integrals. 
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4 ■ 3. Regularization 

The second method does not make use of the formal expression (2.18), but utihzes the 
original expression (2.16) with the regulating parameters —L_ = L_|_ = L. The integral 
involved is, for the j = state. 



rL /-tanh^L 



J-L Jo 

= 2F(i, 1 + iujo; |; tanh^L) tanhL 

= BU, -iuo) - ' F{1, i - iuo; 1 - iuo] cosh'^L) tanhL , (4.12) 

where F{a, b; c; x) is the hypergeometric function, and the last step follows from its 
transformation properties. This then gives 

{<Po\4>o) = -li^o) - ^KiL) , (4.13) 

UJq 

where 

K{L) =D\- iF(l, i - iujQ] 1 - iuoq] cosh'^L) tanhL . (4.14) 

Now it is easily shown, using even the crudest approximation = iujq, that K{L) = 
0(e'^^), whereas (t>l{L) = 0{e^). Thus the second term in ( |4.13| ) (which is guaranteed 
to be independent of L) is shown to be zero when evaluated at L ^ oo. The norm 
of the j = 1 state can also be recovered in this manner, though one needs a better 
approximation for in this case; a way to obtain these better approximations is given 
below. 

The matrix elements are likewise regulated, provided we know the logarithmic 
derivatives D±{uj) when the wavefunction is integrated from x ±oo. These are 
not available for a general potential. Another way to regulate the matrix elements is to 
integrate along a contour in complex x plane, e.g., along the path x = ue^^, u real and 
9 fixed. The matrix element becomes 

/oo 
due''ct>l{ue'')Vo{ue''). (4.15) 
-oo 

It is easily seen that (pKue^^) decays exponentially for sufficiently large 6 and hence 
the integral along the rotated contour converges |T^. As an example, we compute the 
matrix element numerically for j = 1 state using ( 4.15| ) with 6 = 60°. The result agrees 



with the analytic value given by (4.11). In the case of P-T potential, the logarithmic 
derivatives D±[lj) are available analytically. Hence the matrix element can also be 
evaluated by (2.15), but since this result is highly special, we shall not exhibit it here. 
Instead, we go on to a numerical scheme applicable to all potentials with exponential 
tails. 
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4.4- Numerical evaluation and Born series 

For any system, provided the unperturbed wavefunction 0o(a;) is known, the shifts 
are in principle given by (2.15) and (2.16). These involve (a) finite integrals from 
L_ = —L to Lj^ = L, which can be handled numerically in the usual way; and (b) 
surface contributions involving D±. The latter contain all the information from the tails 
of the potential. The lowest approximation Dj-{uj) ^ ±iuj is in general not accurate 
enough, because it multiplies (pl{L) ~ e'^'^°^. 

A general, yet simple way to obtain a better approximation for |x| > L is to use 
the Born approximation. Here V will stand for any potential; by applying the method 
sketched below, we can find the logarithmic derivative of either Vq or Vq + Vi, and hence 
obtain the quantities D±. If V{x) oc e~"^ and the Born approximation is carried to 
mth order, then the remaining error would go as V{x)^~^^ oc e"*^™"*"^)"^', which will be 
sufficiently accurate for dealing with any unperturbed state with 70 < (m + l)a. 

The Born approximation is particularly easy to implement for a potential that goes 
as an exponential. For simplicity we deal with the tail at x = L only. Let 

V{x) = VoJ2cke-'"'\ (4.16) 

k 

where by convention ai < 02 < ■ ■ ■■ The P-T potential is of this form, where ak = 2k /b, 
and Cfc = (-1)^=+^^. 

By iterating (|2.2| ) in powers of V, one finds 

/o(a;) = iu , 

fiix) = VoY.Okiak-2tu)-'e-'''=\ 

k 

f2{x) = Vo EcfcCfc.(«fe - 2tuj)-\ak' - 2tu)-\ak + a^' - 2tuj)-'e-^'"'+''^'^' . (4.17) 

k,k' 

etc. (Here the subscripts on / denote the order of the Born approximation, not LPT.) 
Higher-order terms can be generated readily by algebraic software. All the sums can be 
terminated at some kmax if only accuracy up to 0(e~^^) is required for some finite /?. 

The poles in (4.17) exist only in the Born approximation, but not in the exact 
solution Nevertheless, (4.17) makes it clear that even exponentially small potentials 

% In deriving the Born approximation, one has in effect first taken Vq — > 0, then secondly considered 
say ak — 2iuj — ^ in the resultant expression (4.17). This order of the limits implies that the result is 
only valid for Vbe""*"^ <C \ak — 2iLu\. On the other hand, the exact solution at the position in question 
would refer to taking the limit ak — 2iuj while keeping Vq finite. In this case, there would be no 
pole. 
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can have a significant effect when — Im uj is large, and this is the reason behind the 
string of QNMs — Im AMuok. {k + \) for a Schwarzschild black hole of mass M [32 . 



We have implemented this scheme for the P-T potential, indeed for any potential 
that can be expressed in the form ( ^.161 ) with at = ka. Again for simplicity we deal 
with the situation only on one side, say for the tail as x — > +oo. For this particular form 
of ak spaced evenly in k, the Klein-Gordon equation can be solved easily by substituting 

oo 

(/.(x) = e^"^' ^ rffcC-'^"^ (4.18) 

fc=0 

into (|2l|). One has 



do = 1 



and 



^(icu — ak)dke 

lk=0 



-akx 




(4.20) 



In this example, we have summed four terms with a = 2/6 = 2, and consequently the 
remaining error in the logarithmic derivative calculated is 0(e~^°^'). Incidentally, this 
method, when evaluated at x = L, gives accurate expression for D^{uj) as needed in 
Section O. 



We have used this method to evaluate both the matrix element (2.15) and the norm 
(2.16), taking = L = 5. (In this example, only the positive half line is needed due to 
the symmetry of the potential.) To be precise, the integral over the finite domain (0, L) 
is evaluated numerically, while the surface term is evaluated by the Born series (|4.20|) . 
The result for ui agrees accurately with the result obtained from the two methods 
sketched earlier. 

However, there is still a numerical problem. Take the norm in the j = 1 state as an 
example. Numerical evaluation gives for the two terms in (2.15) 

integral = 187374.578 + 143350. 152z , 
surface term = - 187374.961 - 143350.431? , 

so that there is a loss of 6 significant digits when the two terms are combined. The cause 
of the problem, as before, is the exponential growth of the wavefunction 0o(^) ~ Ae^^°^ , 
so that the asymptotic L dependence of the two terms are respectively ±(A^/2ic(Jo) e^*'^°^, 
where A = 1/2*'^° for the wavefunction (j)Q{x) normalized as in ( f4.9| ). A related difficulty 
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is that the integrand is large and oscillating, which limits the accuracy of evaluating the 
integral. However, these difficulties are readily remedied if we subtract A'^e^^'^"^ from the 
integrand, and add the corresponding term {A'^/2iuJo) [e^^'^°^ — 1) to the surface term. 
Then in this example one finds 

modified integral = - 22.9370946 + 29.1215523i , 
modified surface term = 22.5541253 - 29.4010213? , 

and there is only a loss of 2 significant digits when the two terms are combined. This 
technique can be further refined by removing subasymptotic terms as well. This method 
does not rely on any property of the P-T potential other than the exponential tails. 

The numerical difficulty associated with the exponential growth of the QNM 
wavefunction is exactly the same as the difficulty in the "shooting" algorithm discussed 
in Section |1.2| . This same difficulty, in different guises, always besets numerical solutions 
of QNM problems. Here we have developed an effective method within the realm of 
perturbation theory — but otherwise applicable to any system with exponential tails — 
to tame the problem. The class of problems with exponential tails is sufficiently wide 
for this method to be of interest, especially since there is a dearth of other effective 
methods. With this numerical technique to handle exponential tails, LPT is completely 
formulated for potentials either without tails, or with exponential tails. 

5. Conclusion 

In this paper we have formulated LPT for QNMs. For systems without tails, the 
formalism is no more complicated than for NMs. In fact, there are several advantages: 
the absence of nodes allows simple application to all states, not just the ground state, 
and the possibility that QNMs may not be complete makes alternative methods (e.g., 
generalization of RSPT) less useful. The explicit form of the ffist-order shift is given, 
as well as the most general form of the second-order shift for an arbitrary perturbation. 
When there is a tail that can be expressed as a sum of exponentials, a method is 
developed, based on the Born series, that reduces the calculation to the evaluation of 
integrals, the exponentially large nature of which can be handled by subtracting off the 
leading asymptotic terms. While this is somewhat involved, it is to be stressed that for 
this case no other methods apply in general, not even brute-force numerical integration, 
on account of the need to extract an exponentially small "wrong" solution. Thus the 
technique is likely to be useful. Indeed, this technique has already been employed to 
deal with model perturbation of a black hole pO[] . 

Finally, the generalized norm plays a central role, and in fact has a significance 
beyond perturbation theory. It emerges naturally in the derivation, where the 
integrating factor in ( |2.8D is (pl and not |0oP- However, it is possible to express this 
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same idea in another way, which is possibly more natural and familiar p3[. The idea 



is to write these open systems in terms of a non-hermitian Hamiltonian |25, 26 1 and 
adopt a bi-orthogonal basis which includes a set of left-eigenfunctions dual to the 
right-eigenfunctions 0. Then our generalized norm (0|0) is exactly the same as {4>,(f)), 
where the latter is the conventional inner product which is conjugate linear in the bra 



and linear in the ket. This development will be reported elsewhere p5 
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Figure captions 

Figure la: The trajectory of the lowest eigenvalue in the complex frequency plane as 
Xq is changed from 0.05 to 0.95, for fixed Vq = 100, b = 1, fi = 10, w = 0.1. The 
first-order and second-order perturbation are indicated by the dashed line and the solid 
fine respectively. The triangle and the star mark the positions of the exact eigenvalues 
for xo — 0.05 and 0.95 respectively. The circles show the positions of the exact results 
for other values of xq, and the unperturbed eigenvalue is denoted by the square. 

Figure lb: Same as Figure la but for 1.05 < Xq < 1.50. The triangle and the star mark 
the positions of the exact eigenvalues for xq — 1.05 and 1.50 respectively. 

Figure 2a: The magnitude of the remaining error for the unperturbed value (solid line) , 
first-order perturbation (dashed line) and second-order perturbation (long broken line) 
versus /i, for fixed Vq — 100, b — 1, w — 0.1 and Xq — 0.3. 

Figure 2b: Same as Figure 2a but with xq — 1.4. 



